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Abstract 

We apply periodic orbit theory to a quantum billiard on a torus with a variable 
number TV of small circular scatterers distributed randomly. Provided these scatterers 
are much smaller than the wave length they may be regarded as sources of diffraction. 
The relevant part of the spectral determinant is due to diffractive periodic orbits only. 
We formulate this diffractive zeta function in terms of a N x N transfer matrix, 
which is transformed to real form. The zeros of this determinant can readily be 
computed. The determinant is shown to reproduce the full density of states for 
generic configurations if N > 4. We study the statistics exhibited by these spectra. 
The numerical results suggest that the spectra tend to GOE statistics as the number 
of scatterers increases for typical members of the ensemble. A peculiar situation 
arises for configurations with four scatterers and kR tuned to kR = yo,i ~ 0.899, 
where the statistics appears to be perfectly Poissonian. 

1 Introduction 

Universal level statistics of classically chaotic systems is an asymptotic (i.e. semiclassical) 
property of a spectrum. For homogeneous systems, such as billiards, it appears in the 
high energy limit of the spectrum. Usually there is a preasymptotic regime where the 
statistics reflects the characteristics of classical or quantum origin specific to the system. 

Many authors consider classical diffusion and inhibition of quantum diffusion due to 
localization as the key to understand nonuniversal features of spectra|TJ. In this paper we 
are going to study the problem from a periodic orbit point of view. Admittedly, periodic 
orbit theories based on the semiclassical trace formula Q has not been very successful for 
studies of spectral statistics. The exception is Berry result on the small r limit of the form 
factor ||. This is derived under the so-called diagonal approximation, an approximation 
of somewhat obscure validity^]. 

To simplify the discussion of periodic orbit theories below we will confine ourselves to 
dispersive billiards consisting of one or several disks, all with radius R, inside a rectangle 
or a torus. We thus limit the number of relevant length scales to essentially one. One 
obvious condition for the trace formula to be applicable is kR 3> 1, where k = y/2mE/h is 
the wave number. Non universal effects are expected to be pronounced for intermediate 
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values of kR, and are thus to be sought in the region where the validity of the trace 
formula is dubious. 

Diffraction plays presumably a much more important role in bound system than in 
open and the Gutzwiller formula has to be amended by sums over diffractive periodic 
orbits ||, ||. However, it is not obvious how much the inclusion of diffractive orbits will 
improve the situation. The geometric theory of diffraction is very successful for open 
systems, where it provides small corrections, but its applicability is questionable in the 
penumbra of disk scattering in the very forward direction J?], [|. The basic problem is 
the interference between neutral and unstable orbits accumulating towards them. Their 
respective saddle points are not enough separated and individual eigenstates can hardly 
be resolved by Berry-Keating technique. This should hardly come as a surprise. The 
failure of predicting individual eigenvalues in the semiclassical limit has been expected 
from the very genesis of the trace formula. The pessimism has partly fallen into oblivion 
since the success of the Berry-Keating scheme || |l0| . 

The harsh moral of this discussion is that it seems not an easy task to pursue periodic 
theories as a mean to study how universality may emerge as kR — > do. 

The situation turns out to be very different if the opposite limit (kR <C 1) is con- 
sidered. The disks can now be considered as sources of s-wave diffraction. The unstable 
orbits are replaced by purely diffractive orbits. 

In ref. we studied the small kR limit for the one disk case. The limit is not only 
much easier to deal with than the opposite, but it is also much richer in behavior. In 
particular the A% subspace exhibit a wide range of level statistics in the diffractive region 
which, due to symmetry effects, extend up to kR ~ 4. For (kR ~ 2.40 and kR — > 0) 
the statistics is Poissonian and for (kR ~ 0.899) it is very close to GOE. It approaches 
GOE properly first in the limit kR —* oo. One of the questions we are going to address 
is whether a GOE statistics can be achieved by keeping kR small, but by increasing the 
number of scatterers, and distribute them randomly over a torus. We will thus enter the 
realm of disordered systems. 

In Random Matrix Theories one studies ensemble averages of abstract models of dis- 
ordered systems whereas in Quantum Chaos one usually studies self-averages of chaotic 
systems. There is a need to study how the predictions of Random Matrix Theories is 
realized in more concrete models of disordered systems. We have chosen study spectral 
statistics for individual members of our disordered ensembles, for various values of the 
parameters (N and kR), although nothing prevented us in principle from studying en- 
semble averages. Any study of non- universality and approach to universality for a single 
chaotic systems (utilizing self-averages of the spectrum) will suffer from finiteness of the 
sample. Some type of non-universal effects, like fine wiggles on the formfactor proposed 
by 13] may not relevant when applied to a single systems 1 14], but this is not the kind 
of effect we will be looking at. 

The basic motivation in this paper is conceptual rather than physical, although our 
studies have some bearing on disordered solids and impurity scattering. The system 
considered in this paper have obvious similarities with antidot lattices [|l5|]. However, 
presentday antidot lattices are modeled by rather smooth potentials and they lie in the 
intermediate region kR ~ 2ir (with R suitably defined) where periodic orbit theories are, 
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least to say, cumbersome. 

The outline is as follows. The eigenvalues are recognized as the zeros of the spectral 
determinant, which we will derive within the geometric theory of diffraction. We will thus 
focus on the diffractive determinant (or zeta function), associated with periodic orbits 
with at least one scattering on a diffractive object. This is formulated in section 2.1. This 
object has poles on the real axis so it has to be resummed to tame the divergence caused 
by these poles. This is done in section 2.2. In section 2.3 we discuss some numerical issues. 
In section || we derive the man level density of zeros of the diffractive determinant. In 
section H we compute spectra numerically and study their statistics. 



2 The diffractive determinant 

2.1 Derivation of the diffractive determinant 

In the geometric theory of diffraction^^, [j| |6| the spectral determinant is split up into a 
product 

A(E) = A (E) ■ A G (E) ■ A D (E) , (1) 

where Aq(E) corresponds to the mean level density, the geometric part A G (E) is the 
Gutzwiller-Voros zeta function, possibly amended with the neutral orbits. We will be 
solely interested in the diffractive determinant Ad(E). It has been derived for an non- 
diffractive system supplied with N small disks in ref pl| . 

We thus assume the presence of N small diffractive objects located at r^, where 
1 < k < N whose diffraction constants d k (E) do not depend on the scattering angle. 
We introduce symbolic dynamics by enumerating the disk from 1 to N. The alphabet is 
now {k; 1 < k < N}. The set £Id is defined as the set of all primitive periodic sequences 
of symbols taken from this alphabet. Note that the transition . . . /cjfcj+i . . . does not 
correspond to one trajectory from r/ Cj to r kt+1 , as is usual in symbolic dynamics but all 
trajectories. To clarify this distinction we prefer to talk about periodic symbol sequences 
rather than periodic orbits. 

The diffractive determinant (or zeta function) is now given by 

A D (E)= J] (l-*p) • (2) 

The weight t p is given by 

t p = Y[d ki G G (r ki _ 1 ,r ki ,E) , (3) 

i=i 

where p = k\k2 ■ ■ ■ k np . and /co = k n . Gc(r,r' ' , E) is the non diffractive(i.e. the Green 
function for the original system, before the diffractive objects have been inserted). It 
can, in the semiclassical limit, be decomposed into a sum over all paths (j) from r to r' 

G G (r,r',E)= £ G%\r,r>,E) . (4) 
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where Go is given by the usual van Vleck-Gutzwiller expressions [§]. 

We now turn to our particular system; a rectangle with sides a and b supplied with 
periodic boundary conditions, with N circular disks at positions r& where 1 < k < n. All 
disks have radius R and their diffractive constants are given by 

d[E ) = -4*4^ . (5) 
H$\kR) 

which applies in the limit kR — > 0. We will use as Go the free flight Greens function 
given by 

G Q (r,r',E) = G (r'-r,E)=- l -H^(k\r-r'\) , (6) 

instead of the of usual van Vleck-Gutzwiller, which is obtained by taking the Debye 
approximation of the Hankel functions Hq~\z). The geometric Green function now reads 

G G (r,E)= J2 G o(r + P ,E) . (7) 

p=(ma,nb) 

where the sum runs over all integer m and n. 

Due to a singularity of the Hankel function, this expression diverges if r\2 ~^ 0. We 
define a regularized geometric Green G G (r,E) function by subtracting this singularity. 
The diagonal Green function from a disk to itself is now 

G G (r = 0,E) = - l -J2 H o 1 \kp) , (8) 



and the off diagonal 



G G (r^0,E) = - l -J2Hi ) 1 \k\p + r\) . (9) 



p 



Due to the multiplicative expression for the weights t p @ the diffractive determinant, 
or zeta function, can be computed from the transfer matrix] 17, 15] 



via 



T ij = d{E)-G G {r j -r i ,E) (10) 



£±d{E) = det(l — T) . (11) 

2.2 Making the sums converge and the determinant real 

The sums @ and (|8|) diverge for real E and we will resort to the Ewald summation 
technique as developed in ref. []2(| in order to control the singularities. This procedure 
transforms the diagonal Green function to 

r rn ja - 1 V ex P (Q[i - g 2 /(2£)]) 1 i 

GG ^ E) ~ab „ 4^ (g) + 4 (12) 

g=2TT(m/a,n/a) 
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where I{x) is defined by the integral 

I(x)= / exp(-2xsinhO<2£ • (13) 

Jloslx/Q) 



These expressions are identical to those in [2C], we just keep a and b as free parameters. 



The off diagonal Green function (||) is after resummation 

r l w\ 1 V- < , exp(Q[l-g 2 /(2£)]) 

G g {t 12 ,E) = — 2^cos(ri2-g) 2 £- g 2 ris/O (14) 

g ° 

p 

The derivation of this expression require only slight generalization of the derivation in 
ref. [20], and we omit it. Note that the off diagonal terms are real. 



To get a real expression for the determinant (|ll]) vi simply extract a factor d{E) from 
each row 

A D (E) = {-d) n det{M) , (15) 

where the matrix elements 

M iJo(fcRy + G G (u^J *-J (16) 

3 \ G G (r 12 ,E) i + j 

are all real. 

The energy dependence enter in the Green functions Gg(t, E) and in the renormalized 
diffraction constant 

d=i^M (17) 

where k = y2E. We will in computations artificially fix d and keep the energy dependence 
only in the Green functions, for reasons to be discussed later. 

2.3 Numerical considerations 

The numerical issue is to compute the Green functions with desired accuracy. Through 
the Ewald summation technique each Green function is split up into one sum over the 
dual lattice and one over the initial lattice E There is however no closed expression 
for the function I(x), introduced in Eq. (|l3|). The asymptotic behavior of the function 
I(x) is 

_ exp(Q - x 2 /Q) x 2 /Q - Q 

[ ] ~ Q + xyq [i [xVQ + Qf-' ' [ ' 

As Berry noted, for sufficiently small Q the sum I can De neglected, as far as the 
diagonal Green function is concerned. This is not so for the off diagonal Green functions. 
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They are still small in absolute terms nevertheless significant. We have chosen to include 
the sum J2 !■, and compute the integral (|l^) by the asymptotic formula when appropriate. 
However, for a small number of terms (a number decreasing with increasing energy) the 
integral has to be evaluated numerically. When kr m i n /Q ^> 1, where r m i n the smallest 
interdisk distance, this is no longer an issue. This suggest a small value of Q. One 
the other hand, a large Q is preferred in the dual lattice sum, so the choice of Q is a 
compromise and can be adjusted according to energy. 

The determinant was derived in the limit of small kR. The first correction will involve 
the factor Y\ (kR) / J\ (kR) so our diffractive determinant should work well whenever kR <C 
1. For the N = 1 case the first correction involved the factor Y^(kR) / J^kR), the 
Bessel functions of order 1-3 are suppressed due to symmetry effects pOj]. The diffractive 
approximation then works for slightly higher kR <C 4. Indeed, we found in (l^] that kR 
can be rather close to the limiting kR = 4 (for N = 1) and presumably rather close to 
kR = 1 in the general case. 



3 Mean level density 

In this section we will focus on the mean density of zeros of the determinant A £> (E) as 
given in eqs (11) and (|i~5|). The density of states of the system in the diffractive limit is 
asymptotically given by Weyls expression 

dw = £ • (19) 

This density does not need to be reproduced by the diffractive determinant Ad(E). It 
will turn out that the average density of zeros of Ad(E) will depend on the number of 
scatterers according to 

dW° S = ^f^ dw ■ (20) 

That is, the full spectral density is achieved first when N > 4. We will study spectra for 
fixed values of the parameter d(kR) for reasons to be discussed in the next section. But 
as the result (p0|) does not depend on the value of d, it will also apply to the (physical) 
case where d = d(kR) = d(V2ER) is allowed to vary with energy E. 

Before actually deriving eq ( f20| ) we will make some general comments. 

The reason why we don't resolve the full density of states for N = 1 and N = 2 has a 
simple explanation in terms of symmetries of the system. The wave functions split up into 
the irreducible representations of the respective group and our leading order determinant 
cannot resolve them all. 

If N = 1 the symmetry is C^. I n a coordinate system with origin at the disk there 
is a reflection symmetry with respect to the x and y axis. We only resolve 1/4 of the 
full spectral density, namely those states with even parity with respect both axis. To 
resolve the other subspaces one would need to take higher order terms in the diffraction 
constant, and make a proper desymmetrization. If N = 2 there is a inversion symmetry 
with respect to the point (ri + T2)/2 and we recover 1/2 of the full spectral symmetry. 
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It is obvious that for high enough N one can avoid this kind of symmetry effects. 
It is also obvious that it is possible to construct configurations with large N having 
a high degree of symmetry, like e.g. a regular lattice. Below we will consider generic 
configurations from some random ensemble for which there is no accidental symmetries 
of any kind. We will also assume the absence of exact degeneracies in the unperturbed 
spectrum. 

Now to the derivation of eq ( p0[) . First we show that the mean density of zeros of A^> 
equals the mean density of poles of the same determinant. To show this can essentially 
repeat the arguments in ref |20[| . This was done for a slightly different determinant, but 
the basic mechanism is the same. The difference between the two integrated spectral 
densities is given by [|(]] . 

N D eTOS {E) - N p D oles {E) = -i(Im log Ad(E + ie)) = (21) 

7T 

1 1 

= (Im log det(l - T)) = (Im tr log(l - T)) 

7T 7T 



^ oo 

— Im > 

7T ^ 

r=l 



(trT r ) 



A term trT r is just a product of Hankel functions and is zero on the mean. 

Next we will compute the density of poles of Ad(E). The poles will be located at the 
poles of the Green functions whose density is dw/A- The problem is to determine their 
multiplicity m^. 

So far, we have 

We now study the behavior of det(M) (with M defined in (|l6|)), close to a pole corre- 
sponding to the quantum numbers m, n > 0, that is 2E is close to g 2 = (2ir) 2 ((m / a) 2 + 
(n/b) 2 ). The matrix elements are then (approximately) given by 

M -A J , 4 cos(27rm(x i - Xj)) cos(2irn(y j - yi )) 

M ii -6 i . j a+ ^ 2E-47r 2 (m 2 /a 2 + n 2 /b 2 ) ' [6) 

where we have summed over the four (dual) lattice points 2ir(±m/a, ±n/6). If we intro- 
duce the notation 

ai = 2irmxi/a 



and 



ft = 2irny t /b ' (24) 



X = -^2E-A7T 2 (m 2 /a 2 + n 2 /b 2 ) , (25) 



the latter measuring the (small) distance to the pole, we get 

J 



detM = (-) 7V det(A<5;j — cos(aj — ai) cos(ft,- — ft)) (26) 
A 
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= A"detM . 
A 

The rank of the matrix M is simply the requested multiplicity m n 



detM = (-) n O(\ 
A 



TV-rank(M) 



) - 1/A 



rank(M) 



This rank is the maximum size a matrix having the structure 

Cdj = cos(£j - &) • cos(rjj - rji) 



(27) 



(28) 



can have with a nonvanishing determinant. To explore this problem we introduce three 
other matrices 

CSjj = cos(^- - &) • sin(r/j - 7?i) 



SCi 



sin(^- - &) • cos^- - rji) 



(29) 



SSij = sm(£j - ■ sm(r)j - rji) 

We further introduce the notation CCj to mean the j'th column of CC, and similarly 
for CSj etc. The idea is now to write column CCj as the following linear combination 



CCj = cos(£j — Cj'-i) cos(r]j — rjj-i)CCj-i 

- cos(Cj - Cj-i) sin(r/j - J]j-i)CS^i 

- sin(£j - Cj-i) cos(r/j - r?j_i)SCj_i 
+ sin(^j - sin (??i - ^-i)SSj_i 



(30) 



and similar relations hold for CSj, SCj and SSj. The result can be conveniently ex- 
pressed in terms of matrices. 

u^Tjiij-i , (31) 



where 



u, 



sc- 



and 



and 



T(C,r?) 



/ cos(£) cos(?y) 
cos(£) sin(r/) 
sin(£) cos(ry) 
V sin (6 sin(r?) 



sin(£) sin(ry) \ 
- sin(£) cos (77) 



— cos(^) sin(ry) — sin(£) cos (77) 
cos(^) cos(?y) — sin(£) sin(r/) 

— sin(£) sin(r/) cos(^) cos(r7) — cos(^) sin(r?) 
sin(^) cos(ry) cos(^) sin(ry) cos(^) cos(?y) J 



(32) 



(33) 



(34) 



Please note that the elements of T are scalars, the elements of u are column vectors. The 
external index j has nothing to do with the internal structure of these objects. From the 
definition we have the relation 



T (0 - £j-i,Vj ~ %'-i) T (£j-i - 0-2, ffc-i - = T(^- - Vj - Vj-2) • (35) 
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The basic idea now is to explore whether it is possible to write the column vector 
CC n as a linear combination of the preceding columns CCj (J < n). We address the 
corresponding problem for CS n , SC n and SS n simultaneously, and write 

u„ = T n u„_i = ^ n _iu„_i + (T n - /x n _iE)u n _i , (36) 

where [i n -\ is a multiplier and E is the unit matrix. We carry on this procedure until 
we arrive at 

n-l 

i=2 

where 

S = (T n T n _i---T 2 )-Mn-i(T n _iT n _ 2 ...T 2 )-...M2(T2) , (38) 

or 

n-l 

s = t(£„ - Ci, vn -m)~Yl ^ T te' " ~ ^) • ( 39 ) 

3=2 

So the first n columns of CC are linearly dependent if and only if we can find multipliers 
such that 

S v = j^l , (40) 

The number of multipliers are n — 2 and the number of equations to fulfill is three. 
So for generic parameters £j and 77, the determinant of CC is zero for n > 5. So mjy = 
rankM = min(4,iV). Together with eq. ( p2| ) the announced result (pQ) follows. 

4 Level statistics 

The disks are distributed randomly over the torus according to a uniform distribution. 
We compute spectra for individual members of this ensemble. We choose the lattice 
constants to be a = 1 and b = 2 1//4 , the spectral statistics of the empty torus is than 
perfectly Poissonian; exact degeneracies are avoided. 

The critical parameter is d{kR). As mentioned in the introduction, statistical studies 
suffer for finiteness of the sample. However if R is sufficiently small a sufficiently large 
sample can be obtained with essentially constant kR. But as we also demand that kR < 1 
this would require too large values of k to be numerically tractable. Instead we artificially 
fix kR and and compute the bottom part of the spectrum. We compute around 600 levels 
for each configuration - computations do get a bit tedious for large N. 

We are interested in two measures on the spectra, the integrated level spacing dis- 
tribution P(s) = Jq p(s')ds' where p(s) is the nearest neighbor spacing distribution. 
Secondly we investigate the two point correlation function of levels 

= (T,S((E - Ei)d + e/2) • 5((E - E 3 )d- e/2) ) E (41) 

ij 

= 5(e) + ( ]T S((E - Ei)d + e/2) ■ S((E - E,)d - e/2) ) E = 5(e) + R(e) , 



9 



where the average is taken for a large number of energies. The correlation functions are 
computed over a gaussian window centered at the middle of the sample spectrum, its 
width is about one sixth of the sample size. The results are then smeared with another 
gaussian having width 0.2. 

Below we will restrict our attention to spectral properties of the cases N = 1 and 
N > 4. The reader should bear in mind that for N = 1 only a quarter of the full density 
of states is resolved and the reported result apply to this single subspace. Superposition 
of all four subspace would result in more Poisson-like statistics. 

The underlying spectrum of the empty torus reveals itself clearly in the spectrum for 
large values of d. If N = 1 and d = ±oo the spectrum is Poissonian for a trivial reason, 
the zeros of the determinant have been pushed towards the poles of the Green function 
corresponding to the spectrum of the empty torus, cf. ref. (2(J. If N > 4 and d — * ±oo 
four zeros will be pushed towards each pole. The corresponding limiting integrated level 
spacing distribution is then 

P( s ) = ^ + l(l-exp(- S /4)) . (42) 

This limiting distribution is plotted in Fig. |l| together with results for two different values 
of J (JV = 7). 

For N = 1 and d = the states are, so to say, repelled by the poles of the Green 
function which result in a spectrum exhibiting level repulsion. The level spacings dis- 
tribution is very close to GOE, see 111 and Fig. |. An exact agreement is not possible 



since the eigenvalues are locked between eigenvalues of the integrable torus. This locking 
is released for high enough N In fact, the two point correlation function R(e) shows a 
clearer deviation from GOE than P(s), see Fig. f|. One of the main questions is whether 
GOE can be approached as ./V — ► oo. 

Increasing N only to N = 4, keeping d = (corresponding to kR = 0.899...) 
yields exactly the opposite result. The level spacing distribution appears to be perfectly 
Poissonian, see Fig. [2|. It is known that Poisson-like distribution arises from independent 
superposition of spectra, so one could think that that the determinant (for some unknown 
reason) factorizes. However, the reported distribution agrees better with the Poissonian 
prediction than with the statistics of four superposed Wigner spectra, see Fig. ||. So the 
statistics do appear to be Poissonian. One could also reply that kR = 0.899 ... is to close 
to kR ~ 1 to be physically relevant. However, we know that for N = 1 the diffractive 
approximation is very good close to kR = 0.899 . . . because a pole blows up the element 
in the KKR matrix that corresponds to the diffractive approximation, cf [O], the same 
thing should happen if N = 4, so we are probably talking about a physical effect. 

In Fig. H] we keep d = constant and increase N further. From now on we restrict 
our investigations to the correlation function -R(e). According to the findings for N = I 
we expect it to be a better indicator of deviations from GOE. 

We find that, indeed, the correlator seems to approach that of GOE, for N = 13 it 
already agrees better with GOE than for N = 1. 

However, as we will see, this result is not restricted too any particular choices of kR. 
Next we are going to consider another series of data. Suppose we are increasing N and at 
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\d\ = 


< \d \ < oo 


(i = oo 


N = 1 


« GOE 




Poisson 


N = 4 


Poisson 




4x Poisson 


4 < N < oo 






4x Poisson 


iV = oo 


GOE 


GOE 


GOE 



Table 1: 



the same time decreasing R in such a way that the fraction of the billiard area occupied 
by disks are kept constant: NRlr = C. The corresponding spectra are then studied in 



the neighborhood of some fixed k. For small values of kR we have[21] 



~ lY (kR) 1 „ ,kR, 



We choose arbitrarily k = 2exp7/y(7 and thus d ~ — log N/4tt, and we are led to study 
the sequence 

log N 

d N = -^L . (44) 

The trend is the same, see Fig. [|. The correlation function approaches that of GOE, but 
the approach is of course much slower. 

The conclusions suggested by these studies are summarized in Table |]. The result 
in the lower right corner applies if the limit is approached according to Eq. ( p4| ) or 
something similar. 



5 Discussion 



The emergence of GOE in the limit of many small scatterers will hardly cause any big 
surprise. The interesting thing is that the result has been studied within the framework 
of periodic orbit theory. Admittedly, the periodic orbit was not used directly, that would 
have led to unbearable slow convergence. We chose the underlying system in such a way 
that an efficient resummation could be performed. 

The similarity between the studied system and (disordered) antidot arrays suggests 
that these can be successfully approached from the the diffractive limit rather than from 
the Gutzwiller limit. Higher order terms in the diffraction constant can be treated as 
perturbations. A natural extension of the approach is to apply electric and magnetic 
fields and study transport properties. 



I am grateful to Gabor Vattay for discussions and private lessons on the geometric 
theory of diffraction. This work is a natural continuation of a joint project. This work has 
been supported by the Swedish Natural Science Research Council (NFR) under contract 
no. F-AA/FU 06420-311. 
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Figure 1: The integrated level spacings distribution for N = 7 for two different values of 
d. 
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Figure 2: The integrated level spacings distribution for d = for N = 1 and iV = 4. 
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Figure 3: The integrated level spacings distribution for d = 0, N = 4 compared with the 
Poissonian spectrum and the result from four superposed Wigner spectra. 
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Figure 4: The correlation function -R(e) for d = for variable number of scatterers, 
compared with the GOE prediction. 
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Figure 5: The correlation function R(e) for d chosen according to eq (f44|), for variable 
number of scatterers, compared with the GOE prediction. 
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